Decreasing critical temperature of gas BEC in spatially periodic potential and 
relevance to experiments treated by Mott-Hubbard model. 
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It is shown that the critical temperature of gas Bose- Einstein condensation decreases in deepening 
periodic potential, in contrast to common regularity in a separate potential well. The physical 
explanation of this phenomenon is given. Characteristic scale of potential energies decaying the 
,.0 ' critical temperature is the quantum recoil energy of periodic potential. 

D | The theory represents an alternative and direct approach to the experimental results (C.Orzel et 

Ph ! al Science 291, 2386 (2001); M.Greiner et al, Nature 415, 39 (2002)) obtained with BEC in optical 

l^-} ' lattices and treated as the phase squeezing or Mott transition processes. 

<+H I. INTRODUCTION 

o 

One of the main characteristics of a BEC is the critical temperature of condensation [1] . As is well known, for 
& ' free ideal gas (in thermodynamic limit) it can be simply determined and is presented by formula 
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where m is particle mass, ks is the Boltzmann constant, £ (x) is the Riemann function, N is the number of gas 
particles in volume V . This formula has clear interpretation that the condensation, as a macroscopic quantum 
phenomenon, starts at the temperature when the de Broglie wavelength A c = yj (2ith 2 / mkBT co ) of particle's thermal 

J> ■ motion becomes same order quantity with the distance among the neighboring particles (V/N) 1 ^ 3 ; it ensures the 
existence of perceptible overlapping of neighboring particle wavefunctions (wave-packets), thus making possible the 
manifestation of boson-inductive nature of particles. 

In real situations, when one however has neither ideal nor free gas, aroused the natural question how much and in 



what direction changes the critical temperature compared with the one given by Eq.(I), conditioned respectively by 
interparticle interactions and by existence of trapping external fields (potentials). Interparticle interactions, having 
crucial importance for dynamic properties of condensate, such as stability, excitations and so on, have, fortunately, 
feeble influence on critical temperature of gas condensation [1], [2] (note, that just this feature allows us to use the ideal 
gas model for critical temperature calculations in the following). The influence of the existence of trapping potential 
has been studied rather frequently and the result about the critical temperature can be qualitatively understood by 
means of Eq.(I) or reasonings connected with its physics. Actually, the existence(deepening) of trapping potential 
directly increases the mean density of gas relative to the free state and thereby increases the critical temperature 
too. Here we want to stress, however, that this reason and consequent growth of critical temperature is relevant 
£j \ to the single- well potential and can not be directly implemented for multi-well potentials, in particular for optical 
lattices [3], which attract a great deal of attention thanks to fascinating eventual applications such as matter- wave 
transport [4] and diffraction [5], quantum logic [6], etc. The reason, that in multi-well cases we should exercise some 
caution is that while in single-well case the deepening of potential leads only to drawing nearer all the atoms in the 
trapping potential, in multi-well potentials the reverse tendency exists too. With the deepening of potential, all 
the particles that have preliminarily been in barrier-type regions, fill the neighbor wells, and thereby increasing the 
atomic concentration in well-type regions, decrease it in barrier-type regions. In optical lattices, in addition, the 
mean atomic density is conserved. 

The existence of two opposite tendencies for the gas density, the thickening in well-type regions and respective 
rarefying in intermediate barrier- type regions (see Fig.l), demands a separate considerations of the question [7], even 
more so it may be expected that the critical temperature of condensation in periodic in space potential must be a 
decreasing one depending on the depth of potential, that is it must be straightly opposite to the behavior shown in 
case of single-well potential. To understand the footing of this assertion one should note that the number of atoms 
in barrier-type region in general is less than in well-type region, and besides, the redistribution of atoms between 
barrier-type and well-type regions can be represented as substraction of some (definite) number of atoms from a 
barrier-type region and their addition into a neighbor (left-hand side or right-hand side, it doesn't matter) well-type 
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region. But substraction of a definite quantity from a smaller number is more essential for itself than the addition 
of the same quantity to the larger (or equal) number! That is, the decreasing of the wave-function overlapping in 
barrier-type regions must be more essential for the total gas condensation, than the overlapping increasing for the 
well-type regions, and as a consequence it would lead to mentioned decreasing of the critical temperature. 

In this paper we will afford a quantitative justification of the mentioned assertion about the decaying behavior of 
condensation critical temperature (and as a consequence, of the number of condensed particles) for both deepening 
and/or expanding periodic in space potentials. The form of periodic potential, which will be utilized to this end, will 
be the biparabolic form [8]. On top of that we shall present a simple and convincing physical explanation, based 
on regularities of zonal structure of energy spectrum, why the depth-dependence of critical temperature in standing 
wave trapping potential must be just the contrary to the case of single-well trapping potential. 

We are coming to the problem as in ordinary statistical mechanics, that is, via the familiar relation between the 
number of particles (atoms), temperature and chemical potential [1]. The potential is assumed varying only in one 
direction, meaning that the atomic gas is free in other two directions. For the energy of interacting degree of freedom 
the repeated-zone scheme is used, where the quasimomcntum plays the same role as the ordinary momentum plays 
for free degrees of freedom. A direct consequence of obtained for critical temperature regularity, the eventual 
diminishing of the number of condensed atoms in the deepening potential, also is elucidated:. In addition we calculate 
the number of atoms in separate energy zones as a function of potential depth in order to enlarge the comparison with 
the single-potential case. Finally we work out the dependance of critical temperature on space-period of potential 
and find it to be decreasing too. 

The paper is organized as follows: in Sec. II we shortly represent periodic in one dimension, so-called biparabolic 
potential and Bloch-state solutions for a particle there. In Sec. Ill and IV we outline the theoretical approach, 
illustrate the results of numerical calculations about the critical temperature of condensation, number of condensed 
particles and population of energy zones. Section V contains some comments and conclusions. 



II. BIPARABOLIC POTENTIAL AND DISPERSION RELATION FOR A PARTICLE THERE 

In order to gain an overall picture, let us start and carry the exposition in terms of interactions potential, not 
concretizing, in general, the mechanisms of interaction. However recalling that at the present far off-resonant 
standing wave laser radiation is used to trap condensed dilute atomic gases, we will also write down the expressions 
of potential depth and period for this special case. All the numerical simulations will be performed, having in view 
just this potential. 

Relevant to problem Schrodinger's stationary equation in dimensionless notations is 

*^P. + \W-U(Z)]*(Z) = 0, (2) 

where Z, W and U(Z) respectively are dimensionless coordinate along the periodic potential, total energy relevant to 
this degree of freedom, and potential energy of the particle (atom) there. They arc related with the physical coordinate 
z, energy E, and potential energy V(z) by means of the following formulas: 

Z = 2nz/l, W = E/Er, U(Z) = V(z)/E r , (3) 

where I is the space period of potential, E r =(2-nh) 2 /2ml 2 is the quantum recoil energy of periodic potential, m is the 
particle (atomic) mass. 

As a model of periodic field we choose the biparabolic form [8], 

U(Z) = \—^U + {-l) m k.{Z - rmrf , (4) 

with (m — l/2)7r ^ Z ^ (m + l/2)7r for any integer m — 0, ±1, ±2, .... Here U is the depth of potential, k = 2U/ir 2 is 
a coefficient assuring the continuity of potential, consisting of series of truncated and turned parabolas, as is shown 
in Fig. 2 by a solid line. The potential energy is counted from the bottom of the periodic wells. 

Note that this potential can also be regarded as a good approximation to a potential created for a two-level atom 
by a standing wave of laser radiation, far off-resonant with optical transition. This relevant standing wave potential 
is depicted in Fig. 2 by a dashed-line, and then one has 

2ir/l = 2k, V(Z) = ±\dE L \ 2 /K{uj-uj ), (5) 
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where k = uj/c is the laser wavenumber, d is the dipole matrix element of optical transition of frequency luq, El is 
the amplitude of laser electric field for both traveling waves, composing the standing wave. 
Bloch-solutions of Eq. (2) with potential (4) results in a following dispersion relation: 



cos [2irP] = 1 + 2G n (W)G 22 (W) = 1 - 2G 12 {W)G 21 {W), 



(6) 



where P is the particle's quasimomentum in units 2nh/l( =2hk in laser standing wave case), and the normalized 
energy W is presented in the second equation of (3). Present in (6) functions Gij(W) = 1,2) are expressed 
via the linearly independent solutions Ui(Z) and u 2 (Z) in the well- type region, U\{Z) and u 2 (Z) in the neighbor 
barrier- type region, and their first derivatives at the connective point Z = ir/2 as 
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u 2 {Z) = Zexp 
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and Ui i2 (Z) are obtained from u\ i2 (Z) respectively by means of substitutions Z 2 
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Here $(..., ..; ...) is the confluent hypergeometric function. Let us note that the functions u\ y2 (Z) 

also are real functions of coordinate Z and physical parameters, although they contain imaginary argument — i(Z — tt) 2 
and complex parameter (3 [8]. 



III. CRITICAL TEMPERATURE AND NUMBER OF CONDENSED ATOMS AS A FUNCTION OF 

DEPTH OF PERIODIC IN SPACE POTENTIAL 

Our analysis of ideal gas critical temperature and number of condensed particles (atoms) will be worked out by 
the familiar relation between the mean number of particles, temperature and chemical potential [1], turning from 
the summation over energy spectrum to integration over the momenta. For the energies corresponding to along the 
potential's periodicity degree of freedom, it is convenient use the repeated-zone scheme of energies where each value 
of energy specifies a single value of quasimomentum (see Fig. 3). Quasimomentum naturally plays the role of ordinary 
momentum. In this representation the mentioned relation can be written in form 
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= N. (11) 



V and N are the gas volume and the mean number of particles. The momentum d~P ± (in the free-motion XOY 
plane) is normalized, as the quasimomentum P, in 2ttH/1 units (= 2hk in laser standing wave case), and the energetic 
quantities Wj_, fi and fc^T are normalized in units E r (= (2fifc) 2 /2m in laser standing wave case) as the energy W. 

— >2 

Note that in normalized form W± — P j_ . Integration over the P- momenta is divided into integration over consequent 
zones i = 1, 2, ... . And finally, integration over the negative values is changed into positive-values too, due to which 
the coefficient 2 has appeared in left-hand side of (11). It is to be noted here that the state density is zero at P = 
(as for a free gas) and the twofold repeating of the point P = in (11) does not matter. 
Performing the integration over (free degrees of freedom), instead of (11) we arrive at 
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^f ffM1 _^ IM , (12) 

*- 1 (i-l)/2 

Because of the zero state-density at P = 0, ~P ±_ = (ground state), the critical temperature T c and the number 
of condensed particles N c are afforded to be determined in full analogy with the free ideal gas, that is to consider 
(12) as a determination of /i for a given N, if T T c , but for lower temperatures T ^ T c to fix the value of 
fi(T < T c ) = [i(T = T c ), and consider (12) as a relation for determination of the number of noncondensed particles 
N nc . For the last case, when the gas is cooled lower than the temperature of condensation, the number of condensed 
atoms is [9] 

N c = N - N nc . (13) 

Equations (12), (6)and (13), the second of which determines a one- valued implicit function W(P) in repeated-zone 
picture, are the basis of the reminder of this paper. 

Let us now present the scheme of calculations, which will be used here. To determine T c we first insert the values of 
interest of N, V and I into (12) and consider it as a relation, which determines the chemical potential \i as a function of 
temperature T for a given dispersion relation W{P). Since the seeking functional dependence /x(T) can be obtained 
only by means of numerical integration in(12) and fi sits in integrand, we preliminarily cast the T(/z) dependence and 
hereupon turn to fi(T) function numerically, using the one- valued behavior dependences. In its turn to obtain the T(p) 
dependence we give values to [i and performing the numerical integration in (12) find the corresponding values of T. 
The stage of integration demands of course to fix the depth of the potential too, to have the W(P)— dependence. The 
latter we get by reversing the dispersion relation (6) numerically, where the repeated-zone picture is very convenient. 

Construction of the T(n) dependence in practice starts from some value of /i, greater than the numerical value of the 
particle energy W m \ n for the chosen potential, and gradually move to smaller values of \x. The possible minimal (having 
solution of (12) ) value of fi is fi = W^in, and corresponding to it value of temperature, as solution of (12), just specifies 
the critical temperature T c . Some graphs of the n(T) function for five values of potential depth U, constructed in 
mentioned above way, are depicted in Fig. 4. In accordance with the described procedure of calculations, the value 
of critical temperature for each graph-conditions is determined by the corresponding left-hand side border point of 
graph. Now let us concentrate attention on the fact that the graphs from right to left correspond to lifting values of 
potential depth, which immediately shows the anticipated in Introduction behavior: the critical temperature of ideal 
gas condensation decreases as the potential depth increases. 

It should be noted specially, that the determination of the border points of graphs (and corresponding values of T c ) 
demands a great caution of numerical integration procedure, since the integrand at /i =<! Wmi n (near-border values) 
rapidly increases near the lower edge of the first energy zone, that is near P = 0, where W(P = 0) = W m i n , /j, — W(P = 
0) -> 0, l-exp[(fj,-W(P = 0))/k B T] < 1, and respectively - In [1 - exp [(fi - W(P = 0))/k B T}} > 1. Fortunately, 
the theoretical value of the border (maximal) value of /j, is determined by exact condition ^ max = W m i n , which can 
be fulfilled as precisely as it is possible within the range of computer capability. The minimal accuracy for the 
(dimcnsionlcss) difference has been chosen 10 -8 . Besides, the procedure of integration in (12) was replaced by 
corresponding summation, dividing the first zone quasimomentum range into 10 3 equal parts, and the first of them, 
which is about the P = 0, already into 10 4 equal parts. As a result, the time consumption of our PC for each point 
on the presented graphs was about 30-40 minutes. 

In Fig. 5 we explicitly present the dcpcndancc of gas critical temperature on the depth of periodic potential. It has 
a monotonically decreasing nature everywhere, rather minor in small-depth range and rapid in intermediate-depth 
range of the potential. In the range of greater depths the rate of decreasing , naturally, slows down. 

The number of condensed atoms , as it directly follows from the mentioned behavior of critical temperature, 
decreases with the growth of potential depth too. In addition, since the critical temperature monotonically tends to 
zero as a function of U, for arbitrary chosen low temperature T there is always some potential depth, for which the 
critical temperature becomes equal to T, and starting from which no more condensation takes place. This behavior 
of the number of condensed atoms is shown in Fig. 6. The point where the graph touches the depth axis, just 
corresponds to the peculiar value of potential depth, for which the critical temperature becomes equal to the chosen 
temperature of the gas (0.17o c in this case). For dipper potentials (U > 13) the critical temperature is lower then 
the gas temperature and because of it no condensation takes place. 

Before proceeding to other results obtained, let us turn to other language of presentation, explaining why the 
number of condensed particles decreases with the increasing depths of the periodic potential, which is the equivalent 
of critical temperature decreasing. This representation is based on the behavior of energy spectrum for increasing 
depths of potential and gives, perhaps, more detailed and transparent explanation, than the previous one, based on 
behavior of particle density. It easily explains the studying regularity also in single-potential case. Indeed, in a 
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single-potential the particle has discrete set of energy levels E n (n = 0, 1, ...) of bounded states (only they are relevant 
to the problem under hand) . The probability of population for each excited level, relative to population of the ground 
energy level, is determined by the Boltzmann factor exp [— En/ksT] in thermally equilibrium state. The energy E n 
of each excited energy level can be also interpreted as the distance from ground energy level, choosing without loss 
of generality the energy of the latter zero. The deepening of a single potential, as is well known increases distances 
between energy levels and therefore, between each of them and ground energy level. As a result, the population 
of all excited energy levels decreases. In addition taking into account that the total probability of level occupation 
(population) remains the same unity, we immediately arrive to the result, that the occupation probability of the 
ground energy level must increase due to this deepening of the potential. But the condensed gas is just the particles 
in ground energy level and therefore, the deepening of a single-potential results in increasing of condensed particles, 
or which is equivalent, in increasing of condensation critical temperature. 

In periodic-potential case the situation is, nevertheless, different. The spectrum of a particle has a zonal structure. 
And while the deepening of potential removes from each other the energy bands, it simultaneously narrows all these 
bands, in particular the first one, which plays the dominant role in phenomenon of condensation. In order to easily 
understand the resultant seeking behavior in this situation, it is convenient to present the first energy band as a 
collection of discrete energy levels ( in case of finite, even arbitrary large, number of repeated potential wells this 
approach is absolutely correct ). To the phenomenon of condensation is relevant, by definition, only the ground 
state. Repeating the above discussion about the behavior of level occupation probabilities for deepening potentials 
we see, that the occupation probabilities increase in all excited levels in the first energy band and decrease in all 
upper energy bands. But since the condensation appears at low temperatures, when the low-energy levels are mainly 
populated, the increase of population in excited first-bound levels plays more important role, than the decrease of 
population in upper energy levels. Therefore the increase of potential depth, at least for low temperatures, must result 
in increment of population in the total excited states, including all the excited levels of the first energy band and 
all excited energy bands. This implies that one has decrement of population in ground energy level, relevant to the 
state of condensation. So, regularity of the ground state population for a scries of deepening periodic potentials is 
opposite to the case of a single potential: the number of condensed particles is decreasing and the critical temperature 
of condensation is respectively decreasing. It is not needless to note, that the presented discussion, based on the 
behavior of energy spectrum and relevant Boltzmann function can help us understand qualitatively the behavior of 
number of condensed atoms, or equivalently, the critical temperature in other potential forms too, particularly in 
double well potential. 

In order to give somewhat auxiliary elucidation to the problem under hand, we have calculated the populations of 
each energy band as a function of depth of periodic potential for some, over-critical, temperature, when not only the 
first energy band is effectively populated. The results are presented in Fig.7(point lines). As is seen, the population 
of the first energy band increases monotonically as the depth of periodic potential increases, while the populations 
of the all excited energy bands, vice versa, decrease. It also gives an opportunity to see, that these regularities for 
energy band populations in periodic-wells potential precisely coincide with the regularities of the level-populations for 
single-well potential. This implies, in particular, that the growing behavior of the critical temperature of condensation 
would have been obtained in periodic potential case, if one had taken as a ground level the first energy zone totally, 
which is however incorrect by definition, even for very deep potentials. 

In envisioned circumstances some interest may also present the total number of particles, trapped by the periodic 
potential. That are the particles, energies of which are less than the depth of potential. Because of a zonal structure 
of energy spectrum, in general we should expect a stepwise dependence of this number on the potential depth. The 
results are illustrated in Fig. 8. At the origin (point a\) with U = 0, of course there are no trapped particles. The 
piece ai&i corresponds to situation when the rising potential captures more and more of the first, narrowing at that, 
energy band, naturally increasing the number of trapped particles. Since the energy bands are distanted from each 
other by energy gaps, the pieces corresponding to the second, third, etc. bands are CI2&2, (I363, etc (the upper range 
is not presented), distanted from each other by intermediate pieces 6102, 6203, etc. All the regularities of the graph, 
such as the rate of increasing, saturation, etc., are in agreement with the presented picture, viewing the distribution 
of quantum states within the energy bands and their thermal population as a function of potential depth. 

IV. CRITICAL TEMPERATURE AS A FUNCTION OF PERIOD OF PERIODIC POTENTIAL 

Another side of interest presents the dependence of critical temperature of condensation on the period of periodic 
in space potential. In case of laser-induced standing wave potentials the variation of the period can be obtained 
varying the angle between two waves, composing the standing wave. Denoting by 9 the half of supplementary angle 
between these two wave vectors (see Fig. 9), we get ksw = 2fccos# and respectively I — Tr/kcos9. 8 — corresponds 
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to countcrpropagating case, where the period is minimal. We have found convenient to present the ^-dependence of 
critical temperature as a ^-dependence. The scheme of calculations is the same with the one described in previous 
section. Numerical calculations had been worked out for two {{a)U = 1 and (b)U = 6 ) depths of potential and 
the results are plotted in Fig. 10. The critical temperature, as is seen decreases as the period of potential increases, 
monotonically approaching zero. 

We have not been able to understand this result on the base of picture about the overlapping of particles' material 
waves (see Introduction), but it is easy to do, using the picture of energy spectrum and thermal population of 
energy levels there. Really, the spatial expansion of periodic potential results in narrowing of all energy bands and 
simultaneously in widening of the energy gaps. But it is the same regularity as it was for the deepening potential, 
and therefore the resultant functional dependance also is the same: the critical temperature decreases as the spatial 
period increases. 

Of course, an analogous regularity exists also for the number of condensed atoms, but we do not see the necessity 
of presenting this concrete form. 

V. CONCLUDING REMARKS 

So, the critical temperature of condensation and respectively the number of condensed particles in a trapping 
periodic potential must be less than it is in free state. 

An important question is, of course, how much "dangerous" are these results for experimental situations, typical for 
condensates trapped in optical standing wave potentials. The critical value of thermal energy of free gas condensation 
in typical range of densities 10 11 — 10 12 cm -3 , as is well known is far subrecoil ( is on the order of one percent or less). 
This implies that such a gas may be trapped, in principle, by means of an optical standing wave potential with the 
depth of the same order of magnitude. Looking at the Figs. 5 and 6 we see that this subrecoil range of potentials lies 
close to the origin of horizontal axis, where the lowering of graphs ( critical temperature or/and number of condensed 
particles) is unessential. Moreover, the characteristic parameters of condensation in practice stay unchanged up to 
depths U/E r w 1, that is in so called shallow potentials (note that our scale of recoil energy E r = (2hk) 2 /2m is four 
times greater than the usually using scale h 2 k 2 /2m). However, in deeper periodic potentials with U/E r w 5 -j- 12 (in 
usual scale w 20 4- 50), which are exploited for strong localization of supcr-coolcd condensed gases within separate 
wells, the critical temperature of condensation falls about one order and the gas can get out from the condensation 
state. 

The comparison with the real experimental situation involves another important parameter too: it is the time, 
especially in cases when the gas is filled into the standing wave potential already being in condensed state (obtained 
in magneto-optical traps). The matter is that the above presented theory pertains to thermally equilibrium gas in 
the given external potential field. The gas, finding itself in the periodic field, preliminarily is not in that state, but 
it would tend to that state in the course of time. The dynamics of approaching to the thermally equilibrium state 
is, in general, complicated and requires a separate and accurate consideration, but on qualitative level it may be 
expected, that two processes would determine the transfer of gas into statistical-mechanical equilibrium state in the 
case under consideration. The first one is the process of interparticle collisions, leading to thermal equilibrium, and 
the second one is the tunneling of particle through the potential barrier between neighboring potential wells, providing 
the establishment of zonal structure for particle energy spectrum. Denoting their characteristic times by r co ; and 
i~tun respectively, it will be expected the ultimate thermally equilibrium distribution (with lower critical temperature 
and less number of condensed particles) should be attained in times t ^ r co ;, T tun - 

Let's note that from the viewpoint of our results the distortion of gas BEC-state in a deepening periodic potential 
and its consequent revival in a shallowing one, observed experimentally in [10,?] as a loss and revival of phase, may 
be accounted simply as a respective transfer to noncondensed state and consequent return into the condensed one, 
as is predicted by the above presented theory for a rather slowly deepening and shallowing periodic potential. All 
the more, that the totally distorting potential depths U/E° = 22 (= 5.5 in our £V-scale) in [10] and 44 (= 11 in 
our E^-scale) m [H], he in the region of our curve of Fig5, where the condensed state must be depressed in practice. 
Plus, that in three-dimensional standing wave case [10] the analogous to Fig5 curve must decay more rapidly than 
our, calculated for one-dimensional standing wave. 

The detailed elusidation of results of these experiments in view of presented here theory will be accomplished and 
presented elswhere. 
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FIG. 1. Scheme of gas spacial distribution for (a) a free state, (b) intermediate depth and (c) great depth of periodic potential. 
The mean concentration of gas stays unchanged with the potential depth evolution. 

FIG 2. Biparabolic form of potential (solid line), used in our calculations. For comparison is also presented the respective 
sinusoid (dashed line) of laser standing wave potential. 

FIG 3. Repeated -zone scheme of single particle energy spectrum, which has been used in our calculations. The picture is 
for the case U = 1 (in recoil energy units). 

FIG 4. The temperature dependance of chemical potential fi(N/V,T) for different depths of potential ((a)U = 0, (6)0.5, 
(c)2.5, (d)7,and (e)12) and for N/V = 10 11 cm -3 . Note that the values of /i(vertical axis) are presented relative to minimal 
energy iJmin of particle in periodic potential. Both axis are normalized by recoil energy E r . 

FIG 5. The behavior of critical temperature T c of ideal gas condensation as a function of depth U in external periodic 
potential. The vertical axis T c is normalized by critical temperature T c o of free gas (with N/V = 10 11 cm -3 ), horizontal axis 
U by recoil energy E r , as it is done along all the paper. The discovery of this behavior is the main result of the presented 
study. 

FIG 6. The periodic potential depth dependence of the number N c of condensed particles(atoms). The gas temperature 
is T = 0.1T c o, that is ten times below the free gas critical temperature T c o of condensation (in case under consideration 
fcsTco = 0.002163774E r , or cquivalcntly T c o — 1.03819 * lO -8 ^ if one is concerned with D2 optical transition in sodium atom 

o 

with A = 5896 A). The vertical axis is normalized by the total number N of particles in the gas. 

FIG 7. The periodic potential depth dependence of the numbers N zone of particles, populating first, second, and third zones 
respectively. The latter are given in units of total number N. The gas temperature is chosen ksT — 1.5 and is far over 
the critical value of condensation. The points represent the numbers, irrelevant to the fact are the particles trapped by the 
periodic potential or not, whereas the dashed line represents the number of trapped particles (that is particles with E < U). 
In deep potential regions, when the zone is entirely within the potential, both quantities coincide. As is seen (points) the 
population of the first zone monotonically increases as potential becomes deeper, while the population of the others decreases, 
in full analogy with the regularity of energy- level populations in single -well potential case. 

FIG 8. The periodic potential depth dependence of the number N tr of particles, trapped by the potential. The deepening 
of the potential includes more and more number of zones, thereby increasing the number Ntr- The stepwise behavior of the 
dependence is conditioned by the zonal structure of energy spectrum (for details see the text). 

FIG 9. The geometry of two laser waves, creating the standing-wave periodic potential. 8 — corresponds to the counter- 
propagating case. Variation of 8 changes the spatial period, but keeps the depth constant. 

FIG 10. The behavior of critical temperature T c of ideal gas condensation as a function of angle 8 (see previous Figure). 
The spatial expansion of periodic potential, as the deepening decreases the value of critical temperature of condensation. 
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